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CN \ and scalar fields in brane/bulk hypermultiplets charged under U(l) in the 5D S 1 /Z2 

orbifold model with generic FI terms. A non-trivial VEV of X generates bulk mass terms 
for U(l) charged fields, and their zero modes have non-trivial profiles. In particular, in 
the SUSY-breaking case, bosonic and fermionic zero modes have Gaussian profiles. Such 
non-trivial profiles are useful to explain hierarchical couplings. A toy model for SUSY 
breaking is studied, and it yields sizable D-term contributions to scalar masses. Because 
the overall magnitude of D-term contributions is the same everywhere in the bulk and 
also on both branes, we have to take into account these contributions and other SUSY- 
$_i ' breaking terms to obtain a realistic description. We also give profiles and mass eigenvalues 
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1 Introduction 



For the last several years, 4D N = 1 supersymmetric models with an anomalous £7(1) gauge 
symmetry and non-vanishing Fayet-Iliopoulos (FI) D-term have been studied intensively, and 
many phenomenologically interesting aspects have been obtained. Such anomalous £7(1) gauge 
symmetries, in general, appear in 4D string models. PQ 12 E] This anomaly can be cancelled by 
the 4D Green-Schwarz mechanism. For example, in weakly coupled 4D heterotic string models, 
the 4D dilaton field S is required to transform non-trivially at the one-loop level as 



under the £7(1) transformation with the transformation parameter Ax, where the Green- 
Schwarz coefficient is set as 5gs = trQx/(487r 2 ). Then, a FI term £ is induced by the vacuum 
expectation value (VEV) of the dilaton field S, 



where Mp is the (4D) Planck scale, and K'(S + S) is the first derivative of the Kahler potential 
of the dilaton field. At the tree level, we have K(S + S) = — \n(S + S), and the VEV of the 
dilaton field provides the 4D gauge coupling g\ — \f (ReS). The FI term is suppressed by the 
one-loop factor unless tiQx is quite large. In type I models, other singlet fields, i.e. moduli 
fields, also contribute to the 4D Green-Schwarz mechanism, and a FI term is induced by their 
VEVs.jS] Thus, the magnitude of the FI term is arbitrary in type I models, while its value in 
heterotic models is one-loop suppressed in comparison with Mp. 

Some fields develop VEVs along D-fiat directions to cancel the FI term, and the anomalous 
£7(1) gauge symmetry is broken around the energy scale £ 1//2 , e.g. just below Mp in heterotic 
models. Hence, this anomalous £7(1) gauge symmetry does not remain at low energy scales, 
although its discrete subsymmetry or the corresponding global symmetry may remain even 
at low energy scales. However, an anomalous £7(1) with non-vanishing FI term would create 
phenomenologically interesting effects. Actually, several applications of anomalous £7(1) have 
been addressed. One interesting application is the generation of hierarchical Yukawa couplings 
by the Froggatt-Nielsen mechanism. [4] We could use higher-dimensional couplings as the origin 
of hierarchical fermion masses and mixing angles, where the ratio ^ l f 2 /Mp plays a role in the 
derivation of the mass ratios and mixing angles like that of the Cabibbo angle. [5] 

Another aspect of anomalous £7(1) is that anomalous £7(1) symmetry breaking induces D- 
term contributions to soft supersymmetry (SUSY) breaking scalar masses, jn] These D-term 
contributions are, in general, proportional to the £7(1) charges, and their overall magnitude 
is on the order of other soft SUSY breaking scalar masses. Thus, these D-term contributions 
become a new source of non-universal sfermion masses. This would be dangerous from the 
viewpoint of flavor changing neutral currents (FCNCs). A type of SUSY breaking and the 
mechanism for its mediation by anomalous £7(1) and the FI term has been proposed. (See 
also Ref. jH].) In this type of scenario, D-term contributions can be much larger than other 
soft scalar masses and gaugino masses. For example, in the case that only the top quark has 
vanishing £7(1) charge, stop and gaugino masses can be 0(100) GeV, and the other sfermion 



S — > S + -5 GS A X , 
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Zs=2 5 gs(K'(S + S))M p , 



(2) 
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masses can be 0(10) TeV. This type of sparticle spectrum can satisfy FCNC constraints as 
well as the naturalness problem, and actually, it corresponds to the decoupling solution for the 
SUSY flavor problem. [S] However, with this type of spectrum, we must be careful in regard to 
two-loop renormalization group effects on stop masses, which could become tachvonic.|l(Jj 

Also, the FI term plays a role in the so-called .D-term inflation scenario. |1H W2\ IT3] Other 
interesting behavior has been studied in models with anomalous U(l) and a nonvanishing FI 
term. 

In addition to these 4D models, which were studied several years ago, recently 5D S l /Z 2 
models with the FI term were studied, and interesting behavior was found .[14*]-|23] The S 1 /Z 2 
extra dimension with radius R has two fixed points, which we denote y = and ttR. In general, 
we can put two independent FI terms on these two fixed points, i.e. £o and In Ref. [Hj , 
SUSY breaking in the bulk is studied using the FI term. In Refs. (13 EH HZj , zero- mode wave- 
function profiles are studied, and it is found that their profiles depend on the U{1) charges. 
Hierarchical Yukawa couplings have also been derived. (THl [T7j In these models, the scalar field 
E in the U(l) vector multiplet (in the terminology of 5D N = 1 SUSY ) plays an important 
role. In particular, bulk field profiles are studied systematically in Ref. ^Bj> which considers 
models in which the sum of the FI coefficients vanishes, i.e. £o + £?r = 0. This corresponds to a 
vanishing FI term in the 4D effective theory. 

The interesting behavior exhibited by 4D models mentioned above can exist for non- 
vanishing FI terms. Therefore, here we extend the analysis of Ref. [XBjj to the case with 
£o + £tt 7^ 0. Actually, the SUSY breaking model of Ref. [Hj and the Yukawa hierarchy model 
of Ref. JH] have FI terms for which £ ^ and £ n = 0. Thus, our analysis includes those cases. 
We systematically study the VEV of £ in the generic case that £ an d £rr are independent. 
Then, we consider bulk field profiles. Our result should be useful for further applications. 

In this paper, we study the VEV of £ with generic values of £ and and we consider 
zero-mode wave-function profiles under the non-trivial VEV of £ in the bulk. We also study 
mass eigenvalues and profiles of higher modes. These systematic analyses are the main purpose 
of this paper. In addition, we consider an application to SUSY breaking and examine the 
implication of D-term contributions to SUSY breaking scalar masses. 

This paper is organized as follows. In ^ we systematically study the question of along 
which direction VEVs are developed in models with generic values of Fl-terms. Also, we 
consider the zero-mode profiles of bulk fields under non-trivial VEV of £ as well as profiles 
and mass eigenvalues of higher modes. In JJJ as an example of an application, we introduce a 
toy model for SUSY breaking and discuss the implication of D-term contributions to SUSY- 
breaking scalar masses. Section 4 is devoted to a conclusion and discussion. In the appendix, 
we investigate profiles and mass eigenvalues of higher modes. 

2 5D model with FI terms 
2.1 Setup 

We consider a 5D SUSY model on M 4 x S 1 /Z 2 . Because of orbifolding, SUSY is broken into 
N — 1 SUSY, in terminology of 4D theory. Our model includes a U(l) vector multiplet and 
charged hypermultiplets in the bulk. The 5D vector multiplet consists of a 5D vector field A M , 
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gaugino fields A±, a real scalar field E and a triplet of auxiliary fields D a . Their Z 2 parities 
are assigned in a way consistent with 4D N = 1 SUSY and U(l) symmetry. |14j For example, 
the Z 2 parity of S field as well as A§ is assigned as odd. A 5D hypermultiplet consists of two 
complex scalar fields, <&±, and their 4D superpartners as well as their auxiliary fields. The fields 
$ ± have Z 2 parities ±, and the U(l) charge of the Z 2 even field $ + is denoted by q, while the 
corresponding Z 2 odd field has U(l) charge —q. Z 2 odd fields should vanish at the fixed points; 
e.g. , 

S(0) = E(-kR) = 0, $_(0) = $_(7ri?) = 0. (3) 
Also the first derivatives of Z 2 even fields along y should vanish at the fixed points; e.g., 

<9 y $+(0) = dyQf+frR) = 0. (4) 

In addition to these bulk fields, our model includes brane fields, i.e. 4D N = 1 chiral multiplets 
on the two fixed points, and their scalar components are written 0/ (J = 0, n), where / denotes 
two fixed points. Their £7(1) charges are denoted by qi. As a first step, we do not consider any 
superpotential on the fixed points. 

We consider the following FI terms on the fixed points: 

t{y) = t^{y) + tJ{y-*R). (5) 

Here we simply assume these FI terms. Our analysis is purely classical. One-loop calculations 
are carried out in Ref. j20]-|2S], and they show that the FI terms are generated at the one-loop 
level. § For one-loop generated FI terms, the sum of FI coefficients £0 + £?r is proportional to the 
sum of U{1) charges. Thus, in the case Co + Ctt 7^ 0, we have an anomaly, although such a case 
corresponds to the 4D models mentioned in the introduction, and we are interested in such a 
case. It is not clear whether or not some singlet fields (brane/bulk moduli fields) transform non- 
trivially under the U(l) transformation to cancel the anomaly and thereby generate FI terms 
like Eqs. (JJ) and (J2J). Moreover, in what follows we discuss non-trivial profiles of bulk fields, 
and such profiles might change one-loop calculations of the FI term from the case with trivial 
profiles. Therefore, let us just assume the above FI terms and study their aspects classically. 

Putting all the above consideration together, the scalar potential terms relevant to our study 
are written [TllTol lTK] 

V = \\-d y Y> + £(y) + gq\§ + \ 2 - gq\§_\ 2 + £ gq I \<p I \ 2 5(y - IR)} 2 

1 /=0,7T 

+ \d y $>+ - gqZ<$>+\ 2 + |<9 y $_ + #g£$_| 2 + 3y|$ + $_| 2 + ---, (6) 

where g is the 5D gauge coupling, which is related to g 4 as g 2 = 2tt Rg\. We have taken a 
vanishing bulk mass, m$ = 0, for &±. The first term corresponds to the D term, and the other 
terms correspond to the F terms. Thus, the .D-flat direction is written 

-d y H + i{y) + gq\^ + \ 2 -gq\^A 2 + £ 9qi\h\ 2 5(y - IR) = 0, (7) 

r=o,7r 

§ In Refs. H3 El it is shown that the one-loop generated FI terms include the second derivatives of 
delta-functions, 6"(y) and S"(y — ttR) as well as delta-function forms as in Eq. ©. Here we do not consider 
such FI terms, but we can extend our analysis to the case with 5"{y) and 5"{y — ttR). 
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and the F-flat direction satisfies 



d y $ ± ^gqE$ ± = 0, (8) 
= 0. (9) 

It is convenient to integrate Eq. (JJJ), and the D-flat condition then becomes 

+ + r R dy(gq\<S> + \ 2 -gq\<5>-\ 2 ) + l- £ gq^i? = 0. (10) 
2 Jo 2 /=0 )7r 

Here we have used the boundary conditions E(0) = E(7nR) = for £(?/). Hereafter, we choose 
the signs of the U(l) charges q and qi such that £ + > 0. We also denote the ratio of £ w to 

£o as rl = 

One of important aspect of the present system is that when E develops a VEV, bulk masses 
of bulk matter fields are generated. When SUSY is not broken, the profiles of the zero modes 
of the Z 2 even fields $+ satisfy the equation 

d y $f _ gq (Z)$f = 0; (11) 

that is, the profiles become 



= $W( )exp 



gq dy'(E(y')) 



(12) 



Fermionic superpartners have the same profiles as the scalar fields (|12p. 

If SUSY is broken and SUSY-breaking scalar masses are induced, the zero-mode profiles 
of the bulk scalar components would change, but the profiles of fermionic components do not 
change from their forms given in Eq. (|12p. Furthermore, if SUSY-breaking scalar masses are 
constant in the bulk, zero-mode profiles of scalar components are not changed, but scalar mass 
eigenvalues are shifted from zero by SUSY-breaking scalar masses. Thus, the VEV of E is 
important. In the following subsections, we study the VEV of E and zero-mode profiles as well 
as higher modes systematically. 

2.2 The VEV of E and the shape of zero modes 

First, we consider the VEV of Z 2 odd bulk fields $_ in the SUSY vacuum. The F-flat condition 
requires the VEV of $_ to satisfy the following equation: 



($_(y)) = ($_(0))exp 



-gq f ' dy '<E(</)> 
Jo 



(13) 



However, the Z 2 odd bulk fields should satisfy the boundary condition $-(0) = 0. This implies 
that ($_(?/)) vanishes everywhere, unless (E) has a singularity. Thus, we consider vacua with 
($_(?/)) = in the unbroken SUSY case. Also, for the broken SUSY case studied in §2.2.1, 
we first consider the direction along ($_(?/)) = 0, and then we discuss whether or not such 
a direction induces tachyonic masses for the first excited mode of Z 2 odd bulk fields. Note 
that the zero modes of Z 2 odd bulk fields are projected out. The first excited mode has a 
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KK mass (e.g., of 0(1/ R)). If the absolute value of negative D-term contributions to SUSY- 
breaking scalar (masses) 2 is larger than such KK (mass) 2 , Z 2 odd bulk fields could also develop 
VEVs. However, we show that such first excited modes of $_ have non-tachyonic masses in the 
following section. Thus, the vacuum with ($_(y)) = 0, which is studied in §2.2.1, corresponds 
to a local minimum. 



2.2.1 The case with (E) ^ ( v g, q z > 0) 

First we consider the case that only E develops a VEV. This situation can be realized when 
q > and qj > for the entire bulk and the brane fields except for Z 2 odd bulk fields. 

Varying E, we obtain the minimization condition, <9 2 E — d y £(y) = 0. Its solution is obtained 

as 

d y {E) = £(y) + C. (14) 

Since E is a Z 2 odd field, it should satisfy the boundary conditions E(0) = E(7ri2) = 0. The 
integral constant C is obtained from the following equation: 

:R dyd y ^) = ^-^ + C^R. (15) 



o 

The left-hand side of this equation should vanish, because of the boundary condition of E. 
Thus we obtain 

C = -k±k. (16) 
2nR v ; 

The VEV (E) itself is written 

(E(y)) = ^£ sgn(y) + Cy + ^(sgnfo - vri?) + 1) (17) 

in the region < y < tcR, where the function sgn(y) is defined as sgn(y) = —1,0,1 for 
y < 0, y = and y > 0, respectively. The case with ^ ^ and £ n = is discussed in Ref. [Tl] . 
and the case with £ + 6zr — is studied in Ref. [TH] . Thus, the above case includes those of 
these previous models. 

Only with the above VEV of E, we have the vacuum energy 

2 



V, 



Thus, in general, SUSY is broken. However, in the case with £o + ^ = 0, we have C = 0, 
and SUSY is unbroken. It is notable that even in the broken SUSY case, the vacuum energy is 
constant along the y direction. 

With the above VEV (|17jh we consider the profile of bulk fields. As we show in Eq. f!24|) . 
even for the case C ^ 0, SUSY breaking scalar masses are constant along the y direction. Thus, 
using Eq. (|12jl. we can write the profiles of both bosonic and fermionic zero- modes of the Z 2 
even bulk field as 

*+(y) = Aexp[^£oy] for C = 0, (19) 



$+(y) = Aexp 



2nR ^ >,„ y, 




for C ^ (20) 
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Figure 1: Wave-function profile of the zero mode of $+. These are plotted by setting girR^ /2 = 
1 and varying q or = £ n /£o- 



in the region < y < rrR, where A is a suitable normalization factor. Note here, however, that 
scalar masses are shifted from zero by D-term contributions. The profile (|19|) is obtained in 
Refs. [HUHEl) and several aspects are discussed in Ref. [TH]. The profile (|2U|) is a new result. 
Actually this profile is interesting. This form is a Gaussian profile localized at 

y = j^r^R- (2i) 

These profiles are shown in Fig. ^ Note that we are considering the case q > for all Z2 
even bulk fields. This localization point is independent of the U(l) charges; that is, all of the 
zero modes are localized at the same point. For example, if £0 — £td ah of the zero modes are 
localized at the point y = (ttR)/2. In some models the localization point may be outside the 
region [0, ttR], and in such cases, we see only the tail of the Gaussian profile. Even in the limit 
of a very large FI term, i.e. £o;£tt ~~ ¥ °°> the ratio £o/(£o + £tt) can t» e finite. However, the 
width of the Gaussian profile depends on the U(l) charges q. This should be useful, e.g., to 
derive hierarchical couplings using wave-function overlapping. 

If we include a non- vanishing bulk mass 771$, the zero- modes become localized at different 
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points. For example, the term \d y § + — gg£$ + | 2 in the scalar potential (JUJ) changes as 

\d y $+ - gq^ + \ 2 - \d y <f> + + (m 4 - ^S)$ + | 2 , (22) 

when we introduce a no n- vanishing bulk mass m$. This corresponds to a constant shift of (E) 
by m^/gq. Then, the localization points of the Gaussian profiles depend on the values of m$. 
Hereafter, we set m$ = 0, but the following analysis and results can be simply extended to the 
case with non- vanishing m$ in a manner similar to that for the above constant shift of (E) . 

The scalar zero mode has a SUSY-breaking scalar mass term due to the non-vanishing D 
term in bulk 

-gqC\<$> + \\ (23) 
which is constant along the y direction, as stated above. This scalar zero-mode mass is 

m ° = m ^R- (24) 

Note that this D-term contribution is positive for the present case with q > 0. Similarly, the 
brane fields on both branes have the D-term contribution to SUSY-breaking scalar mass 

m D,i = mi 2nR • (25) 

This D-term contribution is also positive in the present case with qi > 0. There is no zero- 
mode with tachyonic mass. Note that the overall magnitudes of D-term contributions to scalar 
masses are universal, up to £7(1) charges in bulk and both branes. This point is considered 
again in 

Next, we discuss profiles and mass eigenvalues of the higher scalar modes $±. From the 
potential (jOJ), we have a $± square term after E develops a VEV (E). It is written 

-££ = ±(Z(y)-d y (Z))gq\$ ± \ 2 + \d y $ ± Tgq(m±\ 2 

= (d 2 t gq^y) - G?<?) 2 <£> 2 ) *± = -€a A, (26) 

where A± = <9 2 =F gq£(y) — (gq) 2 (^) 2 - The eigenvalues of the operator A± correspond to the 
KK spectrum of $±. Note that the term ± (£(y) — d y (E)) gq\Q±\ 2 is the SUSY-breaking scalar 
mass term. Thus, we have to solve the equation 

A±$ ± + A$± = 0. (27) 

In the appendix, its solutions are analyzed. In this section, we give only results. 

The case with £ + ^ = 0, i.e. C = 0, is discussed in Ref. [THj . (See also the appendix.) 
Their (mass) 2 spectrum is obtained for Z2 even and odd fields as 

ml = ^ + (jS^l, (fc = l,2,...) (28) 

For the case with C ^ 0, the KK spectrum y/X^/M up to the 30th excited mode is shown 
in Fig. of the appendix, where M is the representative scale of FI terms. From the figure, in 
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(k) o 

the case = 0, we see that for $± , mj, = behaves as 

A* - (f) 2 {te) 2 ±l} ; (k»RM) 



(29) 

A fc ~ (f) (2k ± I), (k<^RM) 



(k) 

and in the case = 1, we obtain for $± 

A fc ~ (f) 2 {fe) 2 ±l}> (*>AM) 

A fc ~ (f) 2 (2A;±l), (k <C i?M) 



(30) 



with the fixed parameters ^o/VnR = M 2 / (\/2it 2 g 4 q) in all cases. The case k <C i?M is rather 
non-trivial, because A& is linear in k and it does not depend on R, the size of the extra dimension. 
This behavior can be understood as follows. Equation (|75|l is a Hermite differential equation 
(harmonic oscillator) for R — > oo. Thus, for & <C i?M, we have a spectrum like that of a 
harmonic oscillator. On the other hand, when A is large compared with g^o/R and gC, the 
term with d y &± in Eq. (|75|) can be ignored. Thus, we have A ~ (k/R) 2 . We also note that the 
first excited mode (k = 1) of $_ does not become tachyonic, because A^/M 2 > (2k — 1)/tt 2 
from Eqs. ([29)1 and f)30j) . Thus, the above vacuum is a local minimum. 

The SUSY-breaking scalar mass terms are constants along the y direction. Hence, higher 
fermionic modes corresponding to $± have the same profiles as their scalar partners. In 
Eqs. (J2HI) and (J3DJ), the terms ±M 2 /(2n 2 ) and ±M 2 /(vr 2 ) are the contributions from the SUSY 
breaking scalar masses. The corresponding fermionic modes do not have such contributions; 
that is, for example, for both = and 1, we have 

A* - (f) 2 (M7) 2 > (k»RM) 
X k ~ (f}2A;. (k<^RM) 

In the case that C is small compared with £o/R, we can consider the term with C in Eq. ([70)1 
to be a perturbation of the case with C = 0. For example, the model discussed in ^corresponds 
to such a case. At the first order of such a perturbation, (mass) 2 eigenvalues A& are shifted 
from (j23) as 

x fc 2 , (Mo) 2 r ^fog)^ y (fe/^) 2 + (Mo) 2 /4 
fc i? 2 4 e^-R _ l (jfe/i2)2 + 5(^)74 + ' C3 2l 

Note that for the first-order perturbation in C exactly cancels the D-term contribution 



There is no tachyonic mode for The mass eigenvalues of the corresponding fermionic 
modes are obtained by adding ±gqC for $ ± . 

2.2.2 The case with (£), (</>/) ^ ( 3 g/ < 0) 

To this point, we have considered the case q > and qi > for all of the brane fields and Z 2 
even bulk fields. In such case, only E develops a VEV. Here, we consider another case, that 
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in which there is a brane field with charge satisfying qj < 0. In this case, such a brane field 
develops a VEV along the D-flat direction, and £7(1) is broken. Our analysis does not depend 
on the brane on which the field with the charge qi < exists. For concreteness, we choose this 
brane to be the y = brane. The D-flat condition is satisfied with the VEV 

gqo(\<Po\) 2 = -(^ + U- (33) 

Now, it is convenient to define the effective FI term, 

?(y) = &5{v) + 0{v-*R), (34) 

where £q = £o + 5"?o(|0o| 2 ) and £^ = Note that for these effective FI coefficients, we have 
£o + £4 = 0- For this case, the D-flat direction for (£) is obtained by replacing £(?/) by £'(y) 
and setting C = in Eqs.(fHj) and (fTTj) . which yields 

(£(</)> = ^osgn(y) + ^;(sgn(y - ttR) + 1). (35) 

Since we assume no superpotential for O , SUSY is unbroken along this direction. 

The zero-mode profile of the Z 2 even bulk field is obtained by replacing £(y) by £,'(y) in 
Eq. (Unj), yielding 



91 c , 



(36) 



This is effectively the same as the case studied systematically in Ref. |18j . 

Mass eigenvalues and wave-function profiles of higher modes are the same as those in the 
case C = considered in $ 12.2.11 except that £ is replaced by that is, (mass) 2 eigenvalues are 
obtained as 

ml ^ 2 + (jSf, (37 ) 
2.2.3 The case with (£), ($+) ^ ( 3 q < 0) 

If there is a Z 2 even bulk field whose charge satisfies q < 0, it can also develop a VEV along 
the flat direction. Here, we consider the case that S and the bulk field $ + with charge q < 
develop VEVs. In this case, the .D-flat and F-flat conditions are written 

-^<£)+£(y) + </ 9 <|$ + |) 2 = 0, (38) 
ay($ + )-gq(E)($+) = 0. (39) 

We multiply the former and latter equations here by (£) and ($+), respectively, and compare 
them. We thereby obtain 

^((S) 2 -(|$ + |) 2 )=2(S)£(y). (40) 

Integrating this equation and substituting the result into Eq. (|3*9*|l. we find that the following 
form of (|$+|) is the solution: 

(|$ + (y)l)= r ^— HTT- (41) 

cos[ay + b + c sgn(y) + c n sgn(y - nR)\ 
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0.2 
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0.6 



0.8 



(a) ($ + (y)) (rj 



y/nR 



(<My)> 




(b) (<& + (y)> ( ff =l) 



<S(y)} 




c) <E(y)> (rff = 1] 




(d) (S(y)> (? = i) 



Figure 2: y dependence of the VEVs ($ + (y)) and (E(y)). These were obtained by setting 
gnR^ /2 = — 1, a = a + (m = 0) and varying g or = £ n /€o- 



Moreover, substituting this form of (|$+|) into Eq. (j39)l . we obtain 

(E(y)) = tanfay + 6 + c sgn(y) + c 7r sgn(y - 7T.R)]. (42) 

Actually, the solution constituted by (}4~Tj) and (|4*2^) is obtained in Ref. [TH] for the case £ ^ 
and £„. = 0. The flat direction represented by (j41j) and (J4*2*|) is the simple extension of that 
result to the case with generic values of £ and £„.. Along this direction, SUSY is unbroken, but 
U(l) is broken. Here, Eq. (J38j) requires the constants to satisfy 

2a 



z 2 2 a2 

a = g q A 



-ci = £/■ 



(43) 



The boundary conditions at y = and ttR imply further restrictions. Explicitly, the boundary 
condition <9 y $ + (0) = E(0) = requires b = c„-,and the same condition of y = nR requires 

anR + c + c w = mn, (44) 

where m is an integer. We combine the latter equation in ()43|) and Eq. ()44j) to obtain a further 
equation that must be satisfied by the constant a, and its solution is obtained as 



a± 



vrm ± ym 2 7r 2 — 2tx Rgq{^ + £„.) 
2kR 



(45) 
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Recall that only the bulk field with charge q < can develop the above VEV. Thus, the value 
inside the square root must always be positive. 

Moreover we have to require the VEVs ($+) and (E) to be non-singular in the region 
< y < tcR. This implies that the singular point y = (mr + 7r/2 — Co)/a, where n is any integer, 
can not be in the region < y < itR. The VEVs of Eqs. (14 lj) and (|42j) are shown in Fig. |21 

Now let us discuss the zero-mode profile of the bulk field with the charge q 1 , whose sign is 
opposite to that of the charge q, i.e. qq' < 0. Such a profile is obtained by substituting this 
VEV of £ into Eq. (JJ. This gives 

(y) = $'j 0) (0) cos 197 " 1 [ay + b + c sgn(y) + c^sgn(?/ - ttR)} . (46) 

Next, we would like to study profiles and mass eigenvalues of higher modes. Some analyses 
are given in the last part of the appendix. Here, we consider the simple case that argument of 
inside the tangent function in (J4"2|) is sufficiently small that we can approximate Eq. ()42j) as 

(S(y)> « -^^V + |sgn(y) + |(sgn(y - - !)» (47) 

where we have taken m = 0. This VEV is exactly the same as that given in Eq. (|17jl. Note 
that in the present case, SUSY is unbroken. Actually, the VEV of $ + (j41|) is approximated as 

<I<M2/)| 2 > « ( 48 ) 

Z7TXI 

at the same level of approximation as Eq. ()47|) . This cancels the VEV (£). Thus, the mass 
eigenvalues and profiles of higher modes are obtained in the same way as in ^2.2.11 but in this 
case there is no SUSY breaking scalar mass terms. Thus, the bosonic and fermionic modes 
have the same spectrum, e.g. Eq. (pTTf . 



3 SUSY breaking from the boundary 

To this point, we have studied VEVs and zero-mode profiles systematically in the cases with 
generic values of £ and £„.. Our main purpose has thus been fulfilled. In this section, we 
consider implications of our results for phenomenological applications. 

In any case, bulk fields have non-trivial profiles, depending on their £7(1) charges. This is 
useful, e.g., to explain hierarchical couplings, as in Refs. [TH] and [TT]. For such a purpose, in 
Ref. ^H] zero- mode profiles with the VEV (J4"2|) in §2. 2. 31 are used, and in Ref. pJJ the zero-mode 
profiles (flUjl are used. The profile (|2"Uj) is also useful. 

Another result, which is important, is the D-term contribution to the SUSY-breaking scalar 
masses (J2~4"j) and (J2"5j) . which are obtained in §2.2.1. The overall magnitudes of the D-term 
contributions are the same, up to U(l) charges, everywhere in the bulk and branes, because the 
VEV of the D term is constant along the y direction. In the case of §2.2.1, SUSY is broken by 
the FI terms, and its breaking scale is O((g(£o + ^/R) 1 ^ 2 )- In this section, we discuss briefly 
D-term contributions in two other cases, that in which the D-flatness and F-flatness conditions 
are inconsistent, which leads to SUSY breaking, and that in which SUSY is broken dominantly 
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by another mechanism independent of the FI terms. In H'd. II we consider a toy model for SUSY 
breaking, which is a 5D simple extension of the 4D models studied in Ref. [7]. In H'3.21 we 
extend the analysis of the D-term contributions in the 4D models studied in Ref. jHj to our 5D 
case. 



3.1 SUSY breaking only by the brane field 

Our starting point here is almost the same as that in §2.2.21 We consider the case that all the 
bulk fields have positive U(l) charges and all of the brane fields also have positive charges, but 
only one brane field Xo at U — has a negative U(l) charge. Here, we employ a normalization 
such that this field has the U(l) charge q xo = — 1. If this brane field has no superpotential on 
the brane, there exists the SUSY-flat direction studied in §2.2. 2| i.e., 

d y {E) = (£ - 9(\xo\ 2 ))5(y) + £J(y - irR), (\ X o\ 2 ) = (49) 

Along this direction, U(l) is broken, although SUSY is unbroken. 
Here we assume the mass term of xo in the superpotential to be 

W = m X ox' S(y), (50) 

where x'o is the brane field on y = and has U(l) charge q x < o = 1. This type of mass term can 
be generated by dynamics on the brane, as discussed in the case of 4D models in Refs. |2j and 
[E]. In this case, the dynamically generated mass can also have a U(l) charge, and other types 
of U(l) charge assignments for x'o are possible. Alternatively, we could write the mass term in 
terms of xo alone. For the sake of simplicity of our toy model, we use the above mass term and 
charge assignment. 

With this mass term, SUSY is broken. Suppose that the field xo develops a VEV (xo) = v. 
Then the analysis in §2.2.11 implies that the scalar potential V is minimized by the following 
VEV of E: 

^<£> = £(y)-gi?8{ y ) + C, (51) 
(£> = l^o-gv 2 )sgn(y) + Cy+^Usgn(y-7cR) + l), (52) 

where C = — (£o ~ yv 2 + Ctt) / (27ri?) . The 4D scalar potential with the above brane mass term 
(|5U|) is written 

V 4 = m v + — . (53) 

Minimizing this potential, we obtain v as 

2 £o + 6r 27ri? 2 

v = —m . (54) 

9 9 2 

With this VEV, the following SUSY-breaking parameters are obtained: 

2 

77? 

- dy (Z)+Z(y)-9(\xo\ 2 ) = — (55) 
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F x' 



^0 + ^ 2tiR 

mv = m\ —m z . (5b) 

V 9 9 
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Note that this F-term is induced only on the y = brane. 

With this VEV of E, the fermionic zero mode as well as the corresponding bosonic mode 
has the profile 

9^o 



<$> { °\y) = Aexp 




2m? 



(57) 



That is the Gaussian profile. Using a method similar to that of ^2.2.11 profiles of higher modes 
as well as their mass eigenvalues are obtained. 

In this vacuum, SUSY is broken and SUSY breaking terms are induced by non-vanishing D 
and F terms. All of the bulk gaugino masses are induced by 

de 2 ^W a WJ(y), (58) 
where A is the cut-off scale; that is, a zero-mode gaugino mass M\ is induced as 

M\ ~ ^ = (59) 

The bulk field with charge q has a SUSY-breaking scalar mass due to non-vanishing D term 
given by 

/Tig = qm 2 . (60) 

The brane field has the same D-term contribution to the SUSY breaking scalar mass, m 2 = qim 2 
. Here, the VEV of the S field plays a role to mediate SUSY breaking to bulk fields and brane 
fields on the other brane as D-term contributions to SUSY breaking scalar masses. 
The brane field on the y = brane has the contact term 

^%rl0o| 2 , (6i) 

which induces the following contribution to scalar masses 

2 2 

m v 



< ~ -nr- ( 62 ) 



"<fo a 2 

Similarly, the bulk scalar fields with vanishing £7(1) charge have a scalar mass m% ~ m 2 v 2 / (RA 3 ) 
from the y = brane. 

In the case R — 0(1/ A), we have the relation 

m q : m 0o : \M X \ = m : me : me 2 , (63) 

where e = v/A. If R is larger than 0(1/A), we would have much smaller gaugino masses than 
scalar masses. Also, note that there are many Kaluza-Klein modes below the cutoff scale in the 
case RA > O(10). Thus, for such case radiative corrections from Kaluza-Klein modes to gauge 
couplings |21j, Yukawa couplings |2S] and SUSY breaking terms |2E1IZZ1 may be important. As 
stated above, the mass spectrum of higher modes is obtained similarly to that of §2.2.1. Here 
we do not study the details of this procedure because this is not our purpose. 
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We have a D-term contribution to scalar masses, as in 4D theory. An important point is 
that the D term is constant along the y direction in our toy model. Therefore, bulk fields and 
brane fields on both branes have the same effects. The field £ plays a role with regard to this 
point. Its VEV, in particular d y (H), levels (flattens) the localized D term in the bulk. This is 
the reason that the D-term contribution appears everywhere in the bulk. Therefore, D-term 
contributions, in general, generate charge-dependent sfermion masses everywhere in bulk and 
on both branes. 

A bulk field could play the same role as that of the brane field x'qi an d we would obtain 
the same result. In this case we could also assume that the mass term superpotential exists 
only on the y = brane. We could consider more complicated models, e.g., a model with two 
U(l)s and their FI terms. Such cases will be studied elsewhere, with the purpose of interesting 
aspects, e.g., realistic Yukawa couplings and sparticle spectra. 

3.2 D-term contribution in the 5D model 

In the above, we have assumed that the superpotential (fKUj) breaks SUSY with a combination 
effect of the FI terms. Here, we give a comment on other cases without such a specific superpo- 
tential. With the same system, except the superpotential (J5Uj) . we assume that SUSY is softly 
broken dominantly through some mechanism and a SUSY-breaking scalar mass m Xo of xo is 
induced. Through this SUSY breaking, gaugino masses and soft scalar masses of other fields 
can be induced. In this case, the previous analysis holds, if m is replaced by m Xo . The £ field 
has a non-trivial VEV, as given in Eq. (|55jh and the fermionic zero mode of bulk field has a 
non-trivial profile, as given in Eq. ([57]). Furthermore, bulk fields with charge q and brane fields 
with charge qj on both branes have D-term contributions to scalar masses given by 

2 2 2 2 (ca\ 

m q = q m xo i m i = ( n m xo i (64) 

respectively. Again, the VEV of the £ field plays a role to induce D-term contributions to 
SUSY breaking scalar masses for bulk fields and brane fields on the other brane. 

4 Conclusion and discussion 

We have systematically studied VEVs of E and brane/bulk fields in a 5D model with generic 
FI terms. We have considered three cases: i) the case in which non-vanishing FI terms are not 
cancelled by the VEVs of brane/bulk scalar fields, ii) the case in which FI terms are cancelled 
by the VEV of the brane field, and iii) the case in which FI terms are cancelled by the VEV 
of bulk fields. The non-trivial VEV of £ generates bulk mass terms for U(l) charged fields, 
and their zero modes have non-trivial profiles. In the first case, SUSY is broken, and fermionic 
zero modes as well as the corresponding bosonic modes have Gaussian profiles. In the second 
case, bulk fields have exponential profiles. In the third case, bulk fields have profiles of the 
form (|46|). We also studied profiles and mass eigenvalues of higher modes. These are our main 
results. Non-trivial profiles would be useful to explain hierarchical couplings. 

We also studied a toy model for SUSY breaking, and we obtained sizable D-term contribu- 
tions to scalar masses. The overall magnitude of the D-term contributions are same everywhere 
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in bulk and also on both branes. We have to take into account these D-term contributions and 
other SUSY-breaking terms in a realistic model. 

Our analysis is purely classical. Actually, one-loop calculations were carried out to show 
that the FI term is generated at the one-loop level. In the generic case, U(l) is anomalous. For 
quantum consistency, we would need anomaly cancellation by the Green-Schwarz mechanism 
and/or the Chern-Simons term. For the former, brane/bulk singlet (moduli) fields might be 
transformed as Eq. ([TJ, so that the Green-Schwarz mechanism is effective. However, this is not 
clear. For the latter, the SUSY Chern-Simons term includes the term YiF MN -Fmtv-E] (See also 
Refs. j2H] and j22j .) Thus, a non-trivial VEV of £ might change the gauge coupling in the bulk 
and/or brane, although in most cases U(l) is broken. 

Nontrivial profiles of bulk fields due to (E) might generate FI terms at the one-loop level 
that are different from that generated in the case with (E) = 0. For example, in the case with 
(E) = 0, the bulk field contributes both to £ and £ n with the same weight. [221 UH HZj However, 
suppose the extreme case in which a non-trivial VEV (E) causes the profile of bulk field to be 
localized completely on the y = brane. Then, such a bulk field might contribute to the FI 
coefficient only on the y = brane. We will study this type of consistency elsewhere. 

Also, one-loop calculations show that FI terms can include second derivatives of delta func- 
tions, S"(y) and 5"(y — nR) . [221 CHI UZj Such terms drastically change the behavior of profiles 
on the branes. ^H] We can extend our analysis to the case including second derivatives of delta 
functions. 
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Higher modes 

Here we study in detailed profiles and mass eigenvalues of higher modes with the non-trivial 
VEV of E obtained in §2.2.1 and 2.2.3. 

The case with (£) ^ ( y q, qi > 0) 

First we study the case of §2.2.1. From the potential © we have a &± square term after E 
develops the VEV (E). It is written 

-£g = ±(e(2/)-^(E))^|$ ± | 2 + |^$ ±T ^(E)$ ± | 2 

= -&± ($ =F gqt(y) - (^) 2 (S) 2 ) $ ± = -$W$±, (65) 

where A± = <9 2 =F gq£,{y) — (fi"?) 2 (^) 2 . The eigenvalues of the operator A± correspond to the 
KK spectrum of $±. 
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We solve the eigenvalue equation 

A±$ ± + A$± = 0, 

where A is the eigenvalue of the operator A±. Following Ref. [TB| . we redefine $± as 

$ ± = $ ±e =tOT , 

and substitute this into Eq. ([HHjl . This yields the eigenvalue equation for $± 

{d 2 y ± 2 M (E}^ ± gq {d y (E) - + A) $± = 0, 

with the boundary conditions 



y=0,irR 



y=0,-rrR 



0. 



(66) 



(67) 



(6* 



(69) 



The VEVs <9 y (£) and (E) are given by Eqs. (JUJ) and (fT7)l . respectively. In the bulk (0 < 
y < tcR), the eigenvalue equation (jSHj) becomes 



d y $>± + (A ± # g c) $± = o. 



(70) 



In the case £o + £tt = C — 0, Eq. (ffOjl is solved in Ref. [IB]. The (mass) 2 eigenvalues A^ and 
wave functions of higher modes for Z2 even bulk fields are obtained as 



A fc = ^7 + 



k 2 , (gq£ ) 



R 2 



(A:) 



up to normalization, where 



1 t + k - 



(71) 
(72) 

(73) 



Z 2 odd bulk fields have the same (mass) 2 eigenvalues A&, and the corresponding wave functions 
are obtained as 

up to normalization. 

In the case £ + £tt = C ^ 0, the differential equation (J70J) can be simplified as 



dl§±{z) =f z0,$±(z) + i/$±(*) = 0, 



where z = (Cy + 



and v 



Z = Z0,Zt, 



-2gqC ^ 2 

2 = 20, 

where zq = y v/r^§ an d ^ = ~^f\f^;- By parameterizing 

6 = 2^%M 2 , 



=F |. The boundary conditions (p)9"j) become 



(75) 



(76) 



(77) 
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where I = (0, 7r), is a dimensionless parameter, and M is a representative scale for the FI 
terms, we obtain the boundary parameters zq and Z>R CIS 



2 = _ -RM, = -TqZq, 



where = l^/ho = £*/£o- The solution of Eq. (|75j) is given by 

*±{z) = f?F?(z,v)+ffF£(z,v), 

where fi and are arbitrary constants and 

F±(z,v) = F( T is/2; 1/2; ±z 2 /2), 
F±(z, v) = zF{\/2 =f ^/2 ; 3/2 ; ±z 2 /2). 



(7? 



(79) 



JO) 



Here, -F(a; b; c) is Rummer's hypergeometric function. Requiring that /f* 1 and / 2 cannot be 
zero simultaneously in the boundary condition (|76j). the following equations should be satisfied: 



d z F+(z,u 
F{{z,v 



d z F+(z,u) 



z=z 



F 2 -(z,v) 



d z F+(z,u 
F{{z,v 



d z F+(z,v) 



F 2 -(z,v) 



z=z 



52) 



These determine the eigenvalue u, which also gives A. We numerically solved Eq. (j82|) and 
evaluate the eigenvalue for the typical case shown in Fig. El From this figure, in the case 



0, we see that for ml = A& behaves as 



h ~ (f)*{(^) 2 ±§}> (k»RM) 

(k < RM) 



X k ~ (f) (2*±§), 
and in the case = 1, we obtain for <3>±, 

^ ~ (f) 2 {fe) 2 ±l 
A, ~ (f) 2 (2A;±l), 



(83) 



(A; > RM) 
(k < EM) 



^4) 



with the fixed parameter ^q/^/tiR = M 2 / \\f2-n 2 g^q) in all cases. The case with -C RM is 
rather non-trivial because A& is linear in k and it does not depend on R, the size of the extra 
dimension. This behavior can be understood as follows. Equation (|75j) is a Hermite differential 
equation (harmonic oscillator) for R — ► oo. Thus, we have a spectrum like that of a harmonic 
oscillator for k -C -RM. On the other hand, when A is large compared with g^o/R and gC, the 
term with d y $± in Eq. (f?T)|) can be ignored. Thus, we would have A ~ (k/R) 2 . 

From the spectrum, we can also derive the wave function profile of <3>±. This is shown in 
Figs. E]and|ID These figures reveal that the higher modes of $ + ($_) localize toward (away 
from) the brane(s) with non-zero Fl-term(s). By fixing all the parameters except for q in Zq 
and z n , we also find that the localization profile becomes sharp as the charge q increases. 
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The case with (E), ($+) ^ ( 3 q < 0) 

Next, we analyze the case with q < for only one bulk field $ + , as studied in §2.2.3. We 
consider another bulk field, with charge ±q' in such a configuration. From the potential 
(0) we have a square term after E and $ + develop VEVs (E) and ($+). It is written 

= -$2(^ 2 T Wa(E))-(V) 2 (S) 2 )$' ± = -$2A ± $' ± , (85) 

where A± = <9 2 =F 9l'{dy&)) ~ (fl"?') 2 (^) 2 - The eigenvalues of the operator A± correspond to 
the KK spectrum of 

We again solve the eigenvalue equation The VEVs Sy(E) and (E) in A± are derived 

from Eqs. flUHJ), (jUJ) and (021) • Then, defining z = tan (ay + b + c ), in the bulk (0 < y < nR), 
the eigenvalue equation becomes 

(1 + z 2 )d 2 z V ± {z) + 2zd K & ± (z) - (n ± (n ± + 1) - & ± (z) = 0, (86) 

where n±(n± + 1) = q'(q' ± q)/q 2 and v± = A/a 2 ± g'/g + n ± (n ± + 1). Thus we obtain 

#fc(z) = At/^Cw) + B ± Ql%{iz), (87) 

where -P^(x) (Q^(x)) is the first (second) class associated Legendre function, and y4 ± and B± 
are constants. The mass eigenvalue A can be obtained by requiring that A± and B± cannot be 
zero simultaneously in the boundary conditions, as in the previous case. 
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k 



5 10 15 20 25 30 , 5 10 15 20 25 30 



(c) $+, z = z n = RM/V2 (d) z = z n = RM/y/2 

Figure 3: KK spectrum a/A^/M up to the 30th excited mode for the case zq = RM and 
z n = (7g = l/(47r 2 g|g) and = 0) in (a) and (b), and for the case z Q = z n = RM/^/2 
iho = k* = 1 / (47r 2 (7|g) ) in (c) and (d). The discrete eigenvalues are joined by solid curves. 

From (a) and (b), we see that \ k = (^ff ^(^f ± §} {k > RM), and X k = (f ) 2 (2& ± |) 

(A; < i?M). From (c) and (d), we see that X k = (f) 2 |(^) 2 ±l| (k > RM), and A fe 

= (2k ±1) (k <C RM). The case k <C i?M is rather non-trivial, because X k is linear in k, 

and it does not depend on R, the size of extra dimension. This can be understood from the fact 
that Eq. (fTo*)) is a Hermite differential equation (harmonic oscillator) for R — > oo. We also note 
that the first excited mode of does not become tachyonic, because \ k jM 2 > {2k — 1)/tt 2 , 
from (b) and (d). 
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Figure 4: Wave function profile of <&± up to the 3rd excited mode, with several boundary 
parameters z (z n = 0) defined by Eq. (|7Sj) . The zero mode belongs to $+, and it is flat with 
respect to $ + . To obtain the wave function profile of the original field, $±, all $>±(y) depicted 
here should be multiplied by the Gaussian factor (ffilj) (see Eq. (jHTjl ). 
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Figure 5: Wave function profile of $± up to the 3rd excited mode, with several boundary 
parameters z = defined by Eq. (|75|). The zero mode belongs to $ + , and it is flat with 
respect to $ + . To obtain the wave function profile of the original field, $±, all $±(y) depicted 
here should be multiplied by the Gaussian factor (J20J) (see Eq. (|57j)). 
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